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The gap parameter for color superconductivity is expected to be a gauge invariant 
quantity, at least on the appropriate mass shell. Computing the gap to subleading order 
in the QCD coupling constant, g, we show that the prefactor of the exponential in 1/g is 
gauge dependent off the mass shell, and independent of gauge on the mass shell. 

1. INTRODUCTION 

Quantum chromodynamics (QCD) is the fundamental theory of the strong interactions. 
In strongly interacting matter at large density or, equivalently, large quark chemical po- 
tential /z, asymptotic freedom |1[] tells us that the QCD coupling g, evaluated at the scale 
/i, is logarithmically small, g(fi) <C 1. Then single-gluon exchange, which is attractive in 
the color-antitriplet channel 0, is the dominant interaction between two quarks at the 
^ ■ edge of the Fermi sea. By Cooper's theorem ||, any attractive interaction destabilizes 
the Fermi surface and, at sufficiently small temperature T, leads to the condensation of 
Cooper pairs. In QCD, quark-quark Cooper pairs break the color symmetry. In this 
color superconductor, at T = it costs at least 20 o to excite a particle-hole pair, where 
0o is the value of the superconducting gap function at the Fermi surface. This gap </> 
can be computed using a mean-field approximation to a self-consistent gap equation with 
single-gluon exchange 

Schematically, this gap equation can be written in the form || 



9 2 



dq 



e+ ^ q 




+ + ln( ^ )+ae; 



(1) 



Here, 0^ is the gap function for quasiparticles, taken on the quasiparticle mass shell. 
In general, the gap function <f) + (K) depends on four-momentum = (k ,k). The 
quasiparticle mass shell is defined by 



k = ei , et^^^-ky+\<P + k \\ (2) 

where e\ is the quasiparticle excitation energy. The gap function on the quasiparticle 
mass shell depends on k only, 0^ = + (e^,k). At the Fermi surface, 0+ = (f) . The 
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integration variable q in Eq. ([I]) is the kinetic energy of ultrarelativistic particles. The 
prefactors (,,(3, (3', and a are either constants or depend at most logarithmically on g. In 
weak coupling, g <C 1, the solution of Eq. ([]]) at the Fermi surface is [§-0] 

<f>o = 2b(iexp(-^J i~L + 0(g)} . (3) 

This solution implies ln(/z/0o) ~ l/#> such that the various terms in Eq. ([!]) differ by 
powers of g. The first term ~ ( contains two powers of the logarithm ln(/x/0o). One is 
the well-known BCS logarithm || which arises from the integration over q, 

The other one is the logarithm multiplying (, when the momentum q is near the Fermi 
surface and e+ ~ ~ 0o- This logarithm is special to theories with long-range inter- 
actions, like the exchange of almost static magnetic gluons in QCD ||||. Its origin is 
a collinear singularity when integrating over the angle between incoming and outgoing 
quark momenta in the gap equation. Since In (/x/0o) ~ ^/g 2 , the term proportional to 
( is the leading term on the right-hand side of the gap equation (|IJ); together with the 
prefactor g 2 it is of order O((f>o), and thus matches the order of magnitude of the left- 
hand-side of the gap equation. The value of the coefficient ( determines the constant c in 
Eq. (Q). This constant was first computed by Son 0, 

The second and third terms in Eq. (|TJ) contain subleading contributions to the gap 
equation. These are proportional to a single power of the logarithm ln(/i/0 o ) ~ For 
the term proportional to (3, this logarithm is the BCS logarithm from Eq. (^J). In the 
other term, proportional to (3', this logarithm is explicitly present, but the BCS logarithm 
is absent, because the additional factor e+ cancels the one in the denominator of the 
integration measure, dq/e'f. Together with the prefactor g 2 , these terms constitute a 
contribution of order O(g<po) to the right-hand-side of the gap equation, i.e., an 0(g) 
correction to the leading contribution. 

The contribution proportional to (3 arises from the exchange of non-static magnetic and 
static electric gluons ||. Both types of interactions are short-range: they are screened 
on a distance scale m' 1 , where m g is the gluon mass; m 2 = Nfg 2 /j J 2 /(67T 2 ), Nf is the 
number of quark flavors. Consequently, the collinear logarithm characteristic for long- 
range interactions is absent, and one is left with the BCS logarithm. The contribution 
proportional to (3' arises from the quark self-energy ||||. It is parametrically of the same 
order as the term multiplying (3. The coefficients j3 and f3' in Eq. (H) determine the 
constant b in Eq. @. Son was the first to give an estimate for the constant b, 

> = («) 



with a constant b of order 0(1), which could not be determined in the approach of Ref. 
||. In HH the constant b was computed by solving the QCD gap equation including 
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non-static magnetic and static electric gluon exchange, but without taking into account 
the quark self-energy. In other words, all terms ~ f3 in Eq. (ffl) were collected, but the 
term ~ (3' was neglected. The result is 




6 = 256vr 4 — b' , (7) 



with an undetermined constant b' of order 0(1). In @J^], i.e., without effects from the 
quark self-energy, b' Q = 1. In [§,0, it was shown that the quark self-energy gives rise to 
a term ~ (3' in Eq. (|l|). As this is parametrically also of subleading order in the gap 
equation, it modifies the constant b' Q , 

( 7T 2 +4\ 

6 = exp(^ —J -0.177. (8) 

The fourth term in Eq. ([]]) summarizes sub-subleading contributions. These are of order 
O(g 2 (f>o). It was argued in [pl|,[5|,|IO~| that at this order gauge-dependent terms enter the 
mean-field gap equation for the color-superconducting gap parameter. However, the gap 
parameter is in principle an observable quantity, and thus gauge independent. Therefore, 
it was concluded that one has to go beyond the mean-field approach to compute gauge- 



independent sub-subleading contributions to the gap parameter. It was also shown JTT 
that effects from the finite lifetime of quasiparticles in the Fermi sea influence the value of 
0o at this order. In weak coupling, the terms ~ a in Eq. (fj) are suppressed by one power 
of g compared to the subleading terms and therefore constitute an order 0(g) correction 
to the prefactor b, as indicated in Eq. (§). 

In this note, we first present the gap equation for the quasiparticle and quasi-antiparticle 
gap, including the gauge-dependent terms, and review previous arguments || on why the 
gauge dependence enters at sub-subleading order in the gap equation. These arguments 
were actually incorrect, in that they neglected additional powers of the gluon momentum 
in the gap equation. Naively correcting for these powers, one obtains that the gauge depen- 
dence enters already at subleading order, giving rise to an extra prefactor ~ exp(3^c , /2) 
to the gap parameter (H), where £c is the gauge parameter for the gluon propagator in 
a general Coulomb gauge. This result is similar to other claims made in the literature 
0. Finally, we demonstrate that a careful calculation of the gauge- dependent term on 
the correct quasiparticle mass shell shows that the gauge dependence indeed enters only 
beyond subleading order in the gap equation. Consequently, the gauge dependence does 
not affect the 0(1) result for the prefactor of the gap, as was originally claimed. 

Our units are h — c — ks — 1, the metric tensor is g^ u = diag(+, — , — , — ), and we 
work in a general Coulomb gauge, with gauge parameter ^c- We denote 4-vectors with 
capital letters, = (ko,h), k = |k|, k = k/&. 



2. GAUGE-DEPENDENT TERMS IN THE QCD GAP EQUATION 

For the sake of simplicity, let us focus on the condensation of quarks with two massless 
flavors, forming Cooper pairs with total spin zero. For the discussion of the gauge de- 
pendence, terms in the gap equation arising from the quark self-energy |J can be safely 
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neglected. Thus, the QCD gap equation reads, cf. Eq. (29) of 



•t>t(K) 



6 90 



d 3 q 

(2^)3 



A^(K - Q) 



+ 



QUO) 

<?0 2 -(^) 2 



Tr 
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-e\2 



Tr 



(9) 



Here, 4>h(K) is the gap function for condensation of quarks with chirality h and energy e, 
e e q = yV-eg^+I^Q)! 2 , and 



^(k) 



1 + hjs 1 + e7o7 • k 



(10) 



are projectors onto states with chirality h = ± = r,£ and energy e = ±. The 
Matsubara sum J],™ runs over fermionic Matsubara frequencies q = —i{2n + 1)tiT. We 
first perform the Matsubara sum and then consider the limit T — > 0. The gluon propagator 
in a general Coulomb gauge is, cf. Eq. (30) of |J, 



Aoo(P) 
A ,(P) 
A«(P) 



A z (P)+ec 



p 4 



c 



PoP* 

p 4 



(^•-PiPi)At(p)+e 



c ■ 



p 4 



To simplify the notation in the following, let us write the QCD gap equation @ 
form 

where R e h (K) is the right-hand side of (|9]) for £ c = 0, and introducing P = K — Q 



(11) 

in the 
(12) 



xi{K) = - 9 2 ty: 

6 qo 



d 3 q P»P V 



(2tt) 3 p 4 



4>%.{Q) 



Tr 



I 5 6 ^)^ Tr 



g 2 



^(k)7 M ^(q)7^ 
^(k) 7 ^^(q)7' 



(13) 



The traces can be readily evaluated, 

1 + k q 



l-kq 



Pi ~(k- qf 
Pl~(k + qf 



(14) 



Since the gap equations for right- and left-handed gaps decouple, we shall drop the index 
h in the following. We furthermore focus on the gap equation for the quasiparticle gap 
function (j) + (K), i.e., e — + in Eqs. (||) and (p"3f) . In the next section, we argue that quasi- 
antiparticle contributions to the gap equation for the quasiparticle gap are negligible. We 
then proceed to estimate the magnitude of the gauge-dependent terms in Section 4. 
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3. NEGLECTING THE QUASI- ANTIPARTICLE MODES 

It is permissible to neglect the contribution from the quasi-antiparticle excitations ~ 
fa{Q) in the gap equation ([]), as has been done in Eq. ([I]). The reason is that the quasi- 
antiparticle gap fa(K) is suppressed with respect to the quasiparticle gap fa(K) by at 
least one power of g. This can be checked by simply power-counting the contributions to 
the gap equation for fa (K) using the arguments presented in the introduction. Taking the 
gap function on the quasi-antiparticle mass shell, k = ejjT, i.e., fa(K) = fa(e^,k) = fa, 
and performing the Matsubara sum over q (which puts the internal gap functions on the 
mass shell, fa{Q) = fa(e^,q) = 0^) and the integral over the angle between external 
quark three-momentum k and internal quark three-momentum q, this gap equation has 
a similar structure as Eq. (|l|), with the obvious replacements fa — > fa, — > e k . (The 
contribution from quasi-antiparticles to this gap equation, which is still present in Eq. 
(H), can be discarded, as it does not even have a BCS logarithm.) Since e^T — fM + k, 
the collinear logarithm from almost static, magnetic gluon exchange in Eq. ([!]) is only of 
order 0(1). This leaves the BCS logarithm. At the Fermi surface 

(15) 

With Eq. (H) one obtains fa ~ gcfio, which proves our assertion. 

The complete contribution from quasi-antiparticle excitations to the equation for the 
quasiparticle gap is suppressed relative to the quasiparticle contribution even more than 
just by a single power of g. Neglecting a possible angular dependence entering through 
the traces in Eq. @ and performing the Matsubara sum which picks up the poles at 
go = — ±(a* + <z)> the quasi-antiparticle contribution is proportional to 




(16) 



i.e., at the Fermi surface it is suppressed by a factor ~ gcj) /(2n) relative to the quasipar- 
ticle contribution (which is proportional to the last factor in Eq. (filf )). 

4. THE GAP IS GAUGE DEPENDENT OFF THE MASS SHELL ... 

The argument of 0, as to why the gauge-dependent terms do not enter into the gap 
equation at leading and subleading order, starts with neglecting the terms proportional to 
pi in Eq. (|14D . This is because on the mass shell such terms are of order p\ = (k — q ) 2 ~ 
(e^ ± e+) 2 ~ 0q. One is then left with the spatially longitudinal terms, cf. Eq. (134) of 

[1, 
Tr 



VOX) 7 • P^(q) 7 • P] - " 1± ^ ■ (17) 

J 2 p z 



As k and q approach the Fermi surface, this term vanishes, completing the argument that 
gauge-dependent terms do not enter at leading and subleading order in the gap equation. 

This is, however, a very general argument. So let us explicitly compute X + (K) from 
Eq. (|i~3"D with ([[7]), just to make sure. We had made the seemingly innocuous assumption 
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of neglecting terms ~ p^. In doing so, implicitly we will be working off the mass shell, 
and so find that the gap is gauge dependent. In the next section, we show that this goes 
away on the mass shell. 

Under the present approximations, there is no term depending on q except for the 
denominator in ([13]). Thus, one can immediately perform the Matsubara sum over q 
using 



T j2 WQ) 



ql - (4) 2 26+ 



2np 




(18) 

go 1U ) — q |_ V 1 ) . 

where np{x) = (e x + 1) is the Fermi distribution. For T = 0, np(e^ /T) vanishes, as e+ > 
> 0. Accounting for an additional factor 1/p 2 from Eq. ([13]), one then integrates over 
d 3 q. The integrand does not depend on the polar angle; therefore d 3 q = 2n dqq 2 dcos#. 
Here, cos 9 = k • q. Substituting p = |k — q| for cos#, the following integral appears in 
Eq. (0), 

k+q A {k + qf-p 2 2kq k + q 

\k- q \ p A [k — q) 2 \k — q\ 

As this term is multiplied by {k — q) 2 on account of (|TTD, the singularity at k = q is 
rendered harmless. Inserting this into Eq. (|13|), one obtains 



K J 24 7T 2 J H e+ 



q (k — q) 2 ^ k + q 



k 2B -\k- q \ ■ < 20) 

The momentum dependence of the gap function restricts the g-integration to a narrow 
interval around the Fermi surface, fi — 6 < q < fi + 5, (5 € /j |. At the Fermi surface, 
k = /i, and introducing £ = q — \i one obtains 



(21) 



where we have neglected some terms of order (/^<i5//i<l, and where = y £ 2 + |0^"| 2 . 
As one approaches the Fermi surface, £ — > 0, the second term in parentheses vanishes. 
Consequently, the first term is the dominant one, and it gives rise to a BCS logarithm, 

X» = /*d£^ ~ * 2 • ( 22 ) 



Such a term is of the same order as terms ~ (3 in Eq. ([!]), i.e., of subleading order in the 
gap equation ! It thus affects the prefactor b in Eq. (|3|). A careful analysis shows that b is 
multiplied by a factor exp(3£<7/2). In pure Coulomb gauge, £c — 0, this factor is unity, 
and one obtains the previous result for b, Eqs. (|6]) and (0). A similar dependence of b on 
the gauge parameter was also reported for covariant gauges . The appearance of a gauge 
dependence in those calculations is, however, not surprising, as the gap parameter is not 
calculated on the quasiparticle mass shell, but for Euclidean (imaginary) energies. Away 
from the quasiparticle mass shell, the gap is not an observable quantity, and therefore 
may depend on the choice of gauge. 

The above argument is, however, still incorrect, as terms proportional to pi were ne- 
glected. Therefore, as mentioned previously, we are also computing off the mass shell. 
Restoring these terms and putting k on the quasiparticle mass shell, we shall see that 
the gauge dependence enters neither at leading nor subleading order in the gap equation. 
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5. ... AND INDEPENDENT OF GAUGE ON THE MASS SHELL 



We now compute the gap on the correct mass shell. To this end, the seemingly innocuous 
terms ~ in Eq. flUp must not be neglected, since p = k —q . To perform the Matsubara 
sum over these terms, one uses the following trick: 



</<) 



<P + (Q) 
«) 2 



(23) 



The Matsubara sum can be computed in the standard way. Taking T — > 0, the final result 
is 



<P + {Q) 



<1Q 



2el 



(24) 



Putting everything together and taking the external energy ko on the mass shell, ko = ejj~ , 
one obtains for X^ = X + (e^, k) the following equation: 



xi 



48 n 2 k 2 



dq^- 



2kq 



(k - q)' 



In 



k + q 
\k — q\ 



(25) 



We again introduce the variable £ = q — /x and neglect terms ~(//i< 1. At the Fermi 
surface, k — we are then left with 



12 vr 2 



d£ 



ill 

4/i 2 11 e , 



(26) 



As in Eq. fl21|), the logarithm in parentheses is completely innocuous because of the 
prefactor £ 2 , and can be neglected in the following. 

The difference to Eq. fl21|) is the second term in brackets, which arises from the p\ term 
in Eq. (|14|). Expanding this term around the Fermi surface (£ = 0), one realizes that it 
is of order £ 2 , and therefore negligible compared to the first term in brackets. Naively, 
one would now conclude that there is a large BCS logarithm / d£/e^ ~ ln(/x/0o), just like 
in Eq. (|22|), and the gauge- dependent term contributes to subleading order in the gap 
equation. This is, however, incorrect. The BCS logarithm "builds up" as £ approaches 
the Fermi surface at £ = 0. This build-up requires a coefficient of order 0(1) over the 
whole range of integration. 

In contrast to Eq. (^T|), this is not the case here, because the factor 1 in Eq. fl2"T| ) is 
replaced by 1 — (0 O — e ?) 2 /£ 2 i n Eq. (|26|) . To see the effect on the magnitude of the 
gauge-dependent terms in the gap equation, split the integral over £ into two parts, one 
from to n<f)o, with k ^> 1, and one from K(po to 5. In the first integral 



0(1) . 



In the second integral, 
(0o -£{) 2 2 o e f 



(27) 



(28) 
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Inserting this into (Iffi), one obtains the order-of-magnitude estimate 

~ 9 2 0o • (29) 

In the first integral, the integration measure would in principle give rise to a BCS loga- 
rithm, if the upper limit of integration was large, ~ 5, and not small, ~ 0o- I n the second 
integral, the factor in the numerator of ( pq) cancels with the one in the integration 
measure d£/eg and thus prevents the "build-up" of the BCS logarithm. 

Comparing the final result (|29| ) with the discussion in the introduction, we conclude 
that the gauge- dependent term is obviously of sub-subleading order in the gap equation 
(P. In other words, the mean- field gap equation for the color-superconducting gap in 
QCD is gauge independent to leading and subleading order. Consequently, the gauge- 
dependent terms influence the prefactor of the gap in Eq. ([|) only at order 0(g). The 
present note thus confirms by analytical means the numerical results of Ref. [ffOj]. 
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